On the experimental feasibility of continuous-variable optical entanglement distillation 
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Entanglement distillation aims at preparing highly entangled states out of a supply of weakly entangled pairs, 
using local devices and classical communication only. In this note we discuss the experimentally feasible 
schemes for optical continuous-variable entanglement distillation that have been presented in [D.E. Browne, 
J. Eisert, S. Scheel, and M.B. Plenio, Phys. Rev. A 67, 062320 (2003)] and [J. Eisert, D.E. Browne, S. Scheel, 
and M.B. Plenio, Annals of Physics (NY) 311, 431 (2004)]. We emphasize their versatility in particular with 
regards to the detection process and discuss the merits of the two proposed detection schemes, namely photo- 
detection and homodyne detection, in the light of experimental realizations of this idea becoming more and 
more feasible. 

PACS numbers: 03.67.-a, 42.50.-p, 03.67.Lx, 03.65.Ta 



The ability to distribute entanglement over large distances 
is one of the key pre-requisites for many practical implemen- 
tations of quantum communication schemes. Quite spectac- 
ular experimental progress has indeed been made in recent 
years towards reaching this aim. Several functioning medium- 
distance quantum key distribution schemes have been reported 
and successful tests indicating a violation of Bell's inequali- 
ties have been carried out. 

Needless to say, any mechanism leading to losses and de- 
coherence will eventually deteriorate entangled states into 
merely classically correlated quantum states. Such states may 
then, for example, no longer be useful in the sense that the 
generation of a secure classical key cannot be guaranteed. To 
regain the ability to distribute entanglement in the presence of 
noise, some instance of an entanglement distillation scheme or 
quantum repeaters is required. In such entanglement distilla- 
tion schemes 1 1], highly entangled states are extracted from a 
situation where entanglement is present in only a dilute form. 
In practical optical schemes such methods form one of the 
building blocks towards making long-distance quantum com- 
munication possible. 

A number of distillation schemes have been devised for dis- 
crete degrees of freedom of light (in particular polarization 
degrees of freedom^ and in some instances even been exper- 
imentally realized Q, |j, 01 • Considering the photon-number 
or continuous-variable degree of freedom, in turn, offers an 
interesting alternative to the former setting, allowing also in 
principle the realization of event-ready entanglement distilla- 
tion without the need of destructive post-selection or photon 
counters. In this context Gaussian states and operations are of 
particular interest as they are experimentally relatively easily 
accessible. It came as a surprise, yet, that with Gaussian op- 
erations alone, continuous-variable entanglement can not be 
distilled \y,^,^. This refers to Gaussian input states, and 
manipulation with passive and active optical elements, homo- 
dyne detection, and vacuuniprojections. Fortunately, it was 
then demonstrated in Refs. fs^^ that one only needs to leave 
the Gaussian setting in a single step to break this no-go theo- 
rem. From there on entanglement distillation is indeed possi- 
ble, making use of Gaussian operations, in fact using passive 



linear optics and vacuum projections or homodyne measure- 
ments only. 

As the realizations of the experimentally feasible schemes 
of Ref. lH, 01 is coming closer in view of the techniques 
that have been developed in recent years, it seems worth dis- 
cussing and emphasizing the versatility of this approach with- 
out presenting mathematical details. This is the aim of this 
brief note. As has already been discussed in Ref. (al. pho- 
ton detectors (with relatively small detection efficiency) are 
suitable in the iteration of the scheme, as well as homodyning 
techniques. Under no circumstances, highly efficient photon 
counters with photon number resolution are required. 
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FIG. 1 : One step of the procedure. The dotted rectangle represents 
the local detection. 

The iteration. - In the following we discuss the basic it- 
erative scheme for entanglement distillation of general non- 
Gaussian states developed in Refs. |8, 9] (Non-Gaussian here 
means that the states are non-Gaussian in the photon number 
degree of freedom; hence, the Wigner function of the states is 
non-Gaussian lUillTjl)- Two copies of a non-Gaussian two- 
mode state p serve as inputs to the procedure IIQ1I . These 
inputs are mixed locally at a 50/50 beam splitter each, repre- 
sented by /7bs, leading to 



p' = {Ubs ^ Ubs){p p){Ubs ^Bs)^ 



(1) 



Here, the symbol denotes a tensor product with respect to 
the two identical two-mode states. In turn, is meant as the 
tensor product between the two remote parties (see Fig. [ij. 
Then, one of the outputs of each of the beam splitters under- 
goes a Gaussian measurement. The outputs of the other arms 



of the beam splitters are retained, in case of a successful mea- 
surement event, and are then used as the input of the next step 
of the procedure. This can be done in an iteration, and each 
instance will from now on be referred to as one step of the 
procedure, or a "Gaussification step". One such step is pre- 
sented in Fig. [l] In practice, already a single step alone can 
significantly increase the degrees of entanglement, and this is 
the setting that seems indeed realistic to reach with present 
technology. In Refs. 10131, the action of such an iteration on 
all non-Gaussian two-mode states has been studied in detail. 
In particular, (i) a proof of convergence to Gaussian states has 
been delivered for all pure or mixed initial states (compare 
also Fig. ISj, (ii) the increase of the degree of entanglement 
has been studied, (iii) the increase in the degree of squeezing, 
and (iv) necessary and sufficient criteria have been presented 
to decide when an exact purification of the input state is pos- 
sible. In the single mode case, enhancement of squeezing and 
the loss of non-classical features has been discussed. The ac- 
tion of the scheme under one or a few steps has also been 
studied in detail. 

The non-Gaussian states can be thought of as resulting from 
a number of possibilities; in Refs. 1 8, 9], in particular, an idea 
has been explored that employs the use of photon subtraction. 
Very much related steps for the creation of a non-Gaussian 
state have already been demonstrated experimentally 113 [l3 




to count photons, but merely requires photon detectors dis- 
criminating between presence or absence of photons. This 
is the variant in terms of which most of the results in Refs. 
11301 have been stated. In Ref. f9ll, the implications of detec- 
tion inefficiencies are discussed in great detail. Quite surpris- 
ingly, the scheme is robust with respect to low detection ef- 
ficiencies, lower than the ones that are already available with 
present technology. The degree of entanglement, measured in 
terms of the logarithmic negativity En as a measure of entan- 
glement BItIi . after a number of steps of the procedure, as a 
function of the detection efficiency r] (r] = 1 corresponds to 
perfect detectors) is depicted in Fig.[3l The initial two-mode 
state p is taken in this plot to be p = | V^) (V^| with 



|^) = (|0,0) + e|l,l))/(l + £2), 
for e = 0.95. Other examples are discussed in Ref. 
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FIG. 3: The logarithmic negativity as a function of the detector ef- 
ficiency r] after one (dashed line) and after 10 steps (dotted line) of 
the procedure, for the initial state given in eq. (|3j. The solid line 
represents the logarithmic negativity of the initial state prior to the 
implementation of the procedure. 



FIG. 2: Photon detection variant of one of the two necessary mea- 
surements. 

Photon detection variants. - To realize the feasible distil- 
lation scheme of Refs. |»,¥1 measurements have to be carried 
out, retaining the outcomes corresponding to a Gaussian pro- 
jection. For this key element of the procedure a number of 
approaches are possible. In one variant, the Gaussian projec- 
tion is obtained by projecting onto the vacuum, making use of 
photon detectors which can discriminate between the presence 
and absence of photons. Such a device - which we assume 
to be perfectly functioning for simplicity, an assumption that 
can easily be relaxed - can be theoretically described as im- 
plementing a measurement with Kraus operators |0)(0| ("no 
click") and I - |0)(0| = ^^=1 l^)(^l ("dick"), see Fig.E) 
One step of the procedure, in a successful event of a vacuum 
projection, hence amounts to the transformation of two two- 
mode input states p into 

p" = (O|0(O|([/Bs0t/Bs)(p0p)(t/BS0t/Bs)^|O)0|O)/7V (2) 

where A^ is an appropriate normalization. 

This state will have a higher degree of entanglement, and 
may or may not be used as the input of the next step of the 
procedure. Note that this scheme does not require the ability 




FIG. 4: Balanced homodyne detection variant of one of the two nec- 
essary measurements. 

Homodyne detection variants. - As has been discussed in 
Ref. 1 8], homodyne detection is a feasible alternative to photo- 
detection in the implementation of the basic protocol detailed 
above. The advantage of such an approach is that the detection 
efficiencies are higher for such measurement schemes than for 
photon detection. 

Of particular interest here is balanced homodyne detection, 
see Fig.|4| In homodyne detection, the signal field is combined 
via a beam splitter with a reference field, referred to as local 
oscillator. In balanced homodyne detection with a sufficiently 



strong local oscillator field, field quadratures can be measured. 
In the scheme depicted in Fig. |4] having an additional input, 
followed by passive linear optics, and homodyne detection, 
the above measurement can indeed be realized. This setup is 
in fact the familiar setup to measure to Q-function of a single 
mode. Concerning the choice of the passive optics, see also 
the appendix. This has been discussed in the context of entan- 
glement distillation in Ref. |8], but is as such a well-known 
observation on general measurement schemes 11^ ll^ 12011 . 
Refs. iH la 01 discuss the possibility of using homodyne de- 
tection in order to realize projections onto Gaussian states in a 
deterministic manner in the language of co variance matrices, 
i.e., moments of quadrature operators. Note that other ampli- 
fication schemes, even phase insensitive amplification, would 
in principle also be suitable. 

Balanced homodyne detection leads to the implementation 
of a projection on a coherent states \a) or, in other words, a 
measurement with POVM elements 



{|a)(a|/7r : a G C}. 



(4) 



As a consequence, such a balanced homodyne detection leads 
us to replace the expression eq. (O by 

p'' = {a\ {P\{Ubs Ubs){p&p){Ubs ^ U^sYla) \p)/N 

(5) 
for any complex a, f3 and, again, for appropriate normaliza- 
tion A^. In the language of moments, each of these projec- 
tions refers to one onto a single-mode Gaussian state with 
second moments given by 7 = I, and first moments d = 
(Re(a), Im(a)). If one now accepts only those measurements 
that correspond to complex a and P close to the origin - ef- 
fectively introducing a cut-off- it is effectively as if the orig- 
inal vacuum projection has been implemented. That is, one 
would for a given x > only accept outcomes for which 
in phase space |a| < x. For small a, we can approximate 
\a) = 1 0)+a I l) + 0(a^), so higher order contributions would 
be orthogonal (in Hilbert space) to an arbitrarily good approx- 
imation, depending on x. See also Ref. 12111 for an effective re- 
alization of vacuum projections with homodyne detection. So 
for the purposes of entanglement distillation, these approaches 
are equivalent. With no exception, all results concerning the 
convergence in the iteration towards Gaussian states, the in- 
crease in the degrees of entanglement and squeezing are just 
as applicable as in the previous variant, without modification 
|8, 9], as the postselected state is identical to arbitrary approx- 
imation under appropriate filtering. 

There is a trade off between the achieved rate when filtering 
successful outcomes and the quality of the distilled state: a too 
weak filtering has the same effect as having non-unit detection 
efficiencies in the previous variant, see above. The interesting 
feature of this variant is - as has been pointed out in Ref. 01 - 
that the high detection efficiencies of homodyning techniques 
can be made use of. The disadvantage may be a lower rate in 
the full scheme due to filtering. 

Remarks on other homodyning variants. - Finally, it is 
worth mentioning that a direct homodyne detection of one 
of the output modes of the two parties will again lead to the 
same predictions under appropriate filtering. The action of a 



homodyning measurement is up to displacements the same as 
a local squeezing, followed by a projection onto the vacuum 
state, in the idealized limit of infinite squeezing. For any real 
non-zero squeezing parameter s, we have that 



[U^s.S{s)QS{s)]={), 



(6) 



and hence. 



{{)\{S{s) I)/7bs = (0|(I S{sy)U^^{S{s) S{s)). (7) 

This is, for any finite squeezing parameter 5, the situation in- 
cluding squeezing is the same as if the input state had been 
appropriately squeezed, followed by an inverse squeezing of 
the final outputs. Again, if one does a homodyne detection 
and filters with respect to outcomes close to the origin, the 
results stated in Ref. |8, 91 concerning a general convergence 
to Gaussian states can be applied. Also, all statements in an 
increase of entanglement are valid in the same manner. 




FIG. 5: An instructive plot of the scheme used as a single-mode pro- 
cedure: this plot depicts the Wigner function of the single-mode state 
after zero, one, and two steps of the procedure. The initial state is 
taken to be |0) + s|l))/(l + £^), where again s = 0.95 (taken from 
Ref. 1 9]). 

Non-Gaussian steps. - The above described procedure will 
take non-Gaussian input and yield an output that is a state that 
is closer to being Gaussian. In the experimental preparation 
of the non-Gaussian states, there is a large variety of possible 
approaches. The scheme is sufficiently versatile to be appli- 
cable as such to all non-Gaussian states. If one encounters 
a Gaussian input from a lossy Gaussian channel, modeling 
photon loss and thermal noise |221, then photon subtraction- 
type |[i3l|23|l is one possible reasonable step 10, 01, making 



use of photon detectors and using the outcome associated with 
I — |0)(0| = X]^i |^)(^| but any operation yielding a non- 
Gaussian state is acceptable. A fair figure of merit for an as- 
sessment of the quality of the procedure would than be the in- 
crease of the degrees of entanglement or purity of the channel 
output compared to the final output of the scheme. In spe- 
cial instances, one may also suffer mere classical ignorance, 
due to classical displacements in phase space due to phase 
noise. This kind of phase noise is related to the so-called 
classical noise channel 11221] . This also leads to mixed non- 
Gaussian states, in case the classical weight of the mixing is 
non-Gaussian. Again, the above procedure would be applica- 
ble. In cases where the classical information about the random 
displacements could in principle be retrieved, however, it may 
be advantageous to directly correct for such errors compared 
to resorting to entanglement distillation. In any case, the en- 
tanglement after Gaussification will typically be higher than 
the one of the mixed state before Gaussification (but gener- 
ally smaller compared to the state before mixing). 

Notably, there are still many challenges to be overcome in a 
full experimental realization of such an idea: mode matching 
at the two beam splitters is definitely an issue, which points to- 
wards the requirement of realizing the scheme entirely within 
fibers, to avoid coupling losses ll24ll . The scheme is obviously 
also demanding in that a number of squeezers with the capa- 
bility to create large degrees of squeezing is required. Dark 
counts and inefficient detectors can to some extent be harm- 
ful, the former of which may be significantly diminished by 
appropriate temporal gating. All these obstacles constitute a 
challenge to the experimental realization, but do not render it 
prohibitively difficult. After all, all ingredients of this scheme 
have already been experimentally realized. 

To summarize, in this note we have emphasized the^xperi- 
mental feasibility of the scheme presented in Refs. 1 8, 9]. It is 
notably sufficiently versatile to allow for homodyne detection 
techniques as measurement schemes, instead of photon detec- 
tion. This may be advantageous when it comes to exploiting 
higher detection efficiencies, together with appropriate filter- 
ing. This alternative was discussed already in Ref. |0], and 
holds as well for the results of the later Ref. 121 • Yet, with ex- 
perimental implementations becoming increasingly feasible, 
this point seems worth emphasizing. It is the hope that this 
note fosters further experimental work in this direction. 
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Appendix. - In this appendix, we represent the arguments of 
Refs. 1510,131 for realizing the above POVM elements. In one 
case, we consider a bi-partite system, consisting of parts 1 and 
2. Part 2 embodies exactly one mode, part 1 may consist of 
a number of n modes. In this case, we consider a projection 
onto the vacuum |0)(0| in system 2. In the other case, we 
take a tri-partite system. Here, 3 is initially prepared in |0) (0|, 
corresponding to an empty port. Now 2 and 3, each consisting 
of one mode, undergo a certain beam splitter transformation. 
Then they are measured with homodyne detection. 

It is not difficult to see that the resulting transforma- 
tion of system 1 is identical, up to displacements in phase 
space. The beam splitter transformation for this to be 
true is specified by a matrix S acting on the quadratures 
(x3,p3,X2,P2,^i,i,Pi,i,...,3;i,n,Pi,n) of the tri-partite sys- 
tem. This matrix is found to be 



a . . a 
a —a . 
a a 

—a . . a 



hn 



(8) 



where a = 1/a/2. When acting on Gaussian inputs, for ex- 
ample, both schemes will result in a transformation of the co- 
variance matrix 



7i,2 



^2 Ci^2 



of the system 1 and 2, according to iQl 

71=51- 0^2(^2 +l2)-'Ci,2. 



(9) 



(10) 



The covariance matrix, in turn, is defined for a state p centered 
at the origin as the matrix with entries 



(7i,2)i,/c = 2Rar[pOjOk] 



(11) 



where O = (x2,P2,^i,i,£i^,^.,^i,n,Pi,n). For a more de- 
tailed analysis, see Refs. 
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